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This paper discusses the formulation and development of a trajectory reconstruction 
tool for the NASA X— 43A/Hyper— X high speed research vehicle, and its implementation 
for the reconstruction and analysis of flight test data. Extended Kalman filtering techniques 
are employed to reconstruct the trajectory of the vehicle, based upon numerical integration 
of inertial measurement data along with redundant measurements of the vehicle state. The 
equations of motion are formulated in order to include the effects of several systematic 
error sources, whose values may also be estimated by the filtering routines. Additionally, 
smoothing algorithms have been implemented in which the final value of the state (or 
an augmented state that includes other systematic error parameters to be estimated) and 
covariance are propagated back to the initial time to generate the best-estimated trajectory, 
based upon all available data. The methods are applied to the problem of reconstructing 
the trajectory of the Hyper-X vehicle from data obtained during the Mach 10 test flight, 
which occurred on November 16th 2004. 


I. Introduction 

T he NASA Hyper-X Program is a joint Langley Research Center and Dryden Flight Research Center 
effort to develop a hypersonic airbreathing research vehicle, powered by an airframe-integrated supersonic 
combustion ramjet (scramjet) propulsion system. Major objectives of the program are to acquire data on the 
operation of the scramjet propulsion system and to demonstrate controlled hypersonic flight technologies. A 
secondary goal is to evaluate the experimental and computational techniques used in the design and analysis 
of the Hyper-X vehicle in order to facilitate future hypersonic and launch vehicle programs.^’ ^ 

The flight experiment mission profile is shown in Fig. 1, and takes place in the Western Aeronautical 
Test Range (WATR) over the Pacific Ocean. The experiment begins with a boost phase, where the Hyper- 
X Launch Vehicle (HXLV) stack (shown in Fig. 2), consisting of the X-43A/Hyper-X Research Vehicle 
(HXRV), shown in Fig. 3 and described in Ref. 3, mated to a modified Pegasus booster as described in 
Ref. 4, is dropped from the NASA B-52 aircraft at an altitude of approximately 40,000 ft and ascends to 
a high altitude. At this point, the HXRV separates from the stack, stabilizes, and the scramjet propulsion 
experiment begins, lasting several seconds. Following the engine test, the vehicle follows a descent trajectory 
that includes various maneuvers that are optimized for the post-flight extraction of stability and control 
derivatives and other aerodynamic properties. Three flights were planned, two at a test condition of Mach 
7 and one at Mach 10. The first flight on June 2, 2001, ended with a failure during the boost phase of the 
trajectory. The second and third tests, on March 27 and November 16, 2004, were both successful. Further 
information on these test flights and preliminary results may be found in Refs. 5-7. 

A critical component of the post-flight data reduction and analysis tasks to determine propulsion sys- 
tem performance and the extraction of various aerodynamic coefficients is the reconstruction of the vehicle 
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Figure 1. Hyper— X Flight Trajectory 



Figure 2. Hyper— X Launch Vehicle Stack Configuration 


trajectory and associated atmospheric flight conditions. The trajectory reconstruction process uses in-flight 
telemetry data combined with a variety of redundant measurements to generate maximum-likelihood esti- 
mates of the vehicle position, velocity, and orientation time histories, along with estimates of the uncertainty 
associated with them. The Hyper-X vehicle has many instruments that gather data to be used for this 
process, such as the Inertial Measurement Unit (IMU), Global Positioning System (GPS), and ground-based 
radar tracking data. 

A tool has been developed for the purpose of reconstructing the trajectory of the Hyper-X vehicle based 
on multiple data sources. The reconstruction approach is a modernization of that given in Refs. 11, 12 and 
uses a recursive, minimum variance Altering algorithm to process accelerometer and gyroscope measurements 
recorded by the IMU during the flight, as well as position and velocity measurements from outside sources, 
such as ground based radar tracking stations, to determine a maximum-likelihood trajectory and user defined 
systematic error sources along with an estimate of the accuracy of these parameters. The data are processed 
through the use of an Extended Kalman Filter (EKF)^^ for a forward and backward pass through the 
data. The results of the two Altering passes are combined using the Eraser-Potter smoothing method^^’^^ 
to determine the Anal trajectory estimate as a function of time. The method also allows the user to specify 
which systematic error sources should be estimated, such as an accelerometer bias, and is general with the 
type of outside tracking measurements available. 

The details of the trajectory reconstruction code and its application to the Mach 7 flight of the Hyper-X 
vehicle have been described in Refs. 8-10. This paper discusses the application of the tool to the data 
acquired during the Mach 10 test flight. A brief review of the reconstruction method is provided first. 
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II. Reconstruction Method 

The reconstruction method is based on the Extended Kalman filter/smoother, which is essentially a 
nonlinear recursive least-squares algorithm. The method makes use of a forward filter and backward filter 
to process data acquired during the flight. The forward Alter is a predictor-corrector approach in which 
the trajectory predictions are formed by direct numerical integration of the onboard accelerometer and 
gyroscope measurements, and trajectory corrections are provided by the processing of Global Positioning 
System position and velocity measurements acquired by the onboard receiver. Smoothing algorithms have 
been implemented in which the Anal estimate of the trajectory and its associated covariance are propagated 
back to the initial time to generate the best-estimated trajectory, based upon all available data. The 
estimation process also develops estimates of various IMU error parameters, such as bias and scale factor 
errors. The Anal product of this process is a planet-relative trajectory along with all uncertainties. The 
speciflcs of the Altering and smoothing equations are discussed as follows. 

A. Filtering and Smoothing Equations 

The current section shall be concerned with the problem of estimating the state of a nonlinear system of 
ordinary diAerential equations 

x(t) = f(x(t),v(t),t) (1) 

where x is the matrix of the system states, f(x,v,t) is a nonlinear function that describes the evolution of 
the system state, and v is a random process with covariance matrix Q given by Q = E(vv^) and mean 
V = E(y), where E is the expectation of the random variable. The state estimation procedure requires 
knowledge of independent observations of some quantity that may be related to the value of the state by the 
nonlinear algebraic equations 

y(t) = g(x(f),w(t),t) (2) 

given at discrete time intervals, where y is the value of the observation, g is the nonlinear relationship 
between the states and the observations, and w is a random measurement noise variable with covariance 
matrix R and mean w. 

The extended Kalman Alter (EKF) is a somewhat ad-hoc modiflcation of the linear Kalman Alter in 
which the linear state transition is replaced by a numerical integration of the full nonlinear dynamic model 
in between measurement updates, before applying the Kalman update at each measurement time. In addi- 
tion, the nonlinear relationship between the state and measurement data is used to calculate the residuals 
rather than the linear approximate relationship. The EKF uses the current estimated state as the reference 
conditions for computing the system dynamics matrices required to calculate the Kalman gain matrix at 
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each time step. The algorithm is given 


Xfe = Xfe + Kfc [yfc - g (xfc,Wfe,tfc)] (3) 

where x^ is the predicted value of the state at time increment k, yk is the value of the measurement data, 
and g (xfc, Wfc, tk) is the predicted value of the measurement data at time ti~, following from the propagation 
of the state estimate from time tk~i to time tk- The gain and covariance equations are given by 

Kk = PkCl{CkPkCl + BkRkT^iy^ (4) 

Pk = (I-KfeCfe)Pfe (5) 

where Pfc is the solution of the matrix differential equation 

P{t) = AP + P^A^ + BQB^ (6) 

P{tk-i) = Pk-i (7) 


at time tk, calculated through numerical integration. The matrices A, B, C and D used in the Kalman gain 
calculations are found by linearizing the system dynamics about the current state estimate, x^, so that 


II 

df 

(8) 

m = 

df 

(9) 

n 

II 

dg 

x=Xfc,w=0,t 

(10) 

II 

Q 

dg 

(11) 


Note that in the formulation of the filtering equations, the state x could include both dynamic states 
as well as constant model parameters to be estimated. Likewise, the input noise v could include random 
inputs to reflect the model uncertainty (process noise), or the uncertainty associated with model parameters 
that are not included in the state x. The uncertain parameters from the latter case are commonly known 
as “consider parameters.” The two noise inputs to the system dynamics both seek to degrade the state 
covariance matrix in order to reflect uncertainty in the model, but are different in the sense that the process 
noise accounts for effects that are unknown whereas the consider parameters account for effects that are 
known but uncertain. A Kalman filter designed to include parametric uncertainty is typically known as a 
Kalman-Schmidt filter. 

The EKF algorithm implemented in this paper may be used for both a forward pass and a backward 
pass through the measurement data. The final estimate from the forward pass through the data is used as 
the initial condition for the backward pass. The final trajectory is then computed as a linear combination of 
the forward and backward passes using the Fraser-Potter method. ^^’^4 optimal smoothed solution at a 
time increment k is given by 



(12) 

X/c = Pfc P X/j, + Pf,^ Xhj, 

(13) 


where P is the smoothed covariance estimate, and x is the smoothed state estimate, P/ and P^ are the 
covariance estimates from the forward and backward filter, respectively, and Xf and x^ are the state estimates 
from the forward and backward filter. By inspecting Fqs. (12-13) it is apparent that the smoother combines 
the forward estimate at time tk, which makes use of all measurement data from the beginning of the data 
set to time tk, with the backward estimate at time tk, which makes use of all data from time tk to the end 
of the interval. In other words, the smoother makes use of all available data at each point of interest. 

B. Dynamic Model 

This section discusses the dynamic model that is used for numerically propagating the state estimates 
between the Kalman filter updates. The equations of motion are based upon Newton’s laws of motion, 
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and are developed in a local, non-inertial coordinate system. The inertial acceleration and angular rate 
measurements recorded during the flight are used in place of detailed aerodynamic and propulsion force and 
moment models. An approximate representation of the complete gravitational potential is used to develop 
an expression for the gravitational force acting upon the vehicle. The following summary begins with the 
definition of various coordinate systems. 

1. Coordinate Systems 

This section provides a brief review of the coordinate systems used to develop the nonlinear equations of 
motion used in the trajectory reconstruction process. Firstly, the Earth Centered Inertial (ECI) frame is 
the frame fixed at the center of the Earth, or any other primary body of interest. The fundamental plane 
of this coordinate system is the equatorial plane of the Earth, with the a;-axis, exi, aligned with the vernal 
equinox, the 0 -axis, ezi, in the direction of the planetary axis of rotation, and y-axis, eyi, chosen in such a 
way as to complete a right-handed coordinate system. This frame is obviously not truly inertial, although 
for short periods of time the inertial assumption is a valid and common approximation (see Ref. 15). 

The Earth Centered Earth Fixed (ECEF) frame is the coordinate system whose origin is located at the 
center of the Earth, similar to the ECI frame, but the coordinate axes are fixed with respect to the planet 
surface. The x-axis, gxfj and the z-axis, Bzfi are aligned with the prime meridian and the axis of planetary 
rotation, respectively. This definition of the ECEF frame is consistent with the World Geodetic System of 
1984 (WGS-84). The ECEF and ECI frames may be related by the simple transformation 


exF 


Byf 

^ZF 


GeCI2ECEF 


GeCI2ECEF 


exi 

Byi 


[ ^ZI 

cos A sin A 0 
— sin A cos A 0 
0 0 1 


(14) 


(15) 


The geocentric {GC) frame is one whose origin is located at the center of mass of the vehicle, and whose 
0 -axis, ezCi is aligned in the direction of the origin of the ECI (and ECEF) frame. The x-axis of the this 
frame, bxcj is perpendicular to the 0 -axis and in the plane defined by §zf and ezci and oriented such that 
the scalar product of exc and e^F is positive. The GC frame may be related to the ECEF frame by the 
transformation 
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(17) 


The geocentric frame is sometimes called the topocentric or local-vertical frame. The ECI, ECEF, and GC 
frames are shown in Fig. 4. 

A frame similar to the GC frame is the geodetic (CD) frame, sometimes known as the topodetic, local- 
horizontal, or North-Earth-Down (NED) frame. This frame is one in which the origin is located at the center 
of mass of the vehicle, and the x-axis, bxd, points in the direction tangent to the surface of the Earth, in the 
direction of the North pole. The 0 -axis, Szd, is aligned perpendicular to the x-axis and oriented downward 
and is perpendicular to the line tangent to the surface of the planet, while the y-axis, byd is shared with 
the GC frame. The CD frame may be related to the GC frame by the transformation 


f exD 


exG 1 

1 ^ 

eye 

> = GgG2GD < 
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( 19 ) 
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Figure 4. ECI, ECEF, and GC Frames 


where e = 4>gd — 4>, with 4>gd being the geodetic latitude of the vehicle’s sub-latitude point, which is the 
projection of the vehicle position to the surface of the Earth along the ezD axis. The geodetic latitude is 
then the angle between the ezo axis and the equatorial plane of the Earth. The GC and GD frames are 
shown in Fig. 5. 

An alternate measure of the location of the sub-latitude point is the geocentric latitude, 4>gc, which is the 
angle from the equatorial plane of the Earth to the line connecting the sub-latitude point and the center of 
the Earth. Fig. 5 shows the differences between the declination, geocentric latitude, and geodetic latitude 
for a planet with exaggerated eccentricity. Note that the three angles are equal for a spherical planet. Fig. 5 
also shows the physical meaning of the geodetic altitude, hgd, and the magnitude of the radius vector, r - 
two variables that will be used presently. 


Table 1. Summary of WGS— 84 Earth Model Constants 


Constant 

Value 

R® 

2.092565 -10^ ft 

e© 

8.1819191 -10-2 

M 

1.4076444 -lO^® ftVs^ 

•h 

1.0826230 TO"® 


4.1780741 TO"® °/s 


The relationship between the geocentric and geodetic latitude is quite simple, and takes the form 

tan (j)gc = (l - e^) tan 4>gd (20) 

The transformation between declination and geodetic latitude, however, is more complicated. Vallado^® 
gives the formulas 


tan 4>gd 


hgd 


tan (j) - 


Ret, el, sin< 


rcoscp^l — el sin^ 
r cos 4> i?0 


^l-el sin^ (j)gd 


( 21 ) 

(22) 
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Figure 5. Declination, Geocentric Latitude, Geodetic Latitude, Geodetic Frame, and Geocentric Frame 


Sofair^^’^® determined an exact solution of Eqs. (21-22), (j)gd{r,(j)) and hgd{r,<j)). The solution is explicit, 
and is free of singularities at the poles and the equator, although the solution space is confined to a region 
outside a spheroid fixed at the center of the Earth with semimajor axis of approximately 43 km. This 
restriction is of no concern for the present problem as the region is well within the surface of the Earth. 
The transformation of position from geodetic to geocentric coordinates is much less complicated, and may 
be determined directly from Eqs. (21) and (22). The WGS-84 values of R<^, C 0 , and other constants that 
describe the shape, mass, and rotation of the Earth are summarized in Table 1. 

Lastly, the body {B) frame is fixed to the center of mass of the vehicle with the Sxb axis directed through 
the nose of the vehicle, the ezB axis oriented perpendicular to the §xb axis, in the plane of symmetry, with 
the positive direction taken to be toward the bottom of the vehicle, and the byb axis completing the right- 
handed system so that it points out the right side of the vehicle. The B frame is shown in Fig. 3. The 
orientation of the B frame is typically expressed by a 3-2-1 Euler angle sequence relating the GD frame to 
the B frame. The transformation is given by 
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where 4', 0, and 4> are the azimuth, elevation, and bank angles, respectively, also commonly referred to as 
the yaw, pitch, and roll angles. 

Although the Euler angle sequence offers an ease of visualization, other parameterizations of the vehicle 
orientation are preferable for implementation in the equations of motion and filtering algorithms. The 
most advantageous of these parameterizations are the Euler parameters, which are defined in terms of the 
Euler angles as^° 
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The use of Euler parameters leads to a representation of the attitude transformation matrix of the form 


GgD2B = 


So + e?-e|-e§ 2(6162 + 6063) 2(6163-6062) 

2(6162-6063) el - el + el -el 2(6061+6263) 
2(6163 + 6062 ) 2 ( 6263 - 6061 ) 


6 n Cl 69 


(29) 


The Euler parameters constitute a four parameter set that describe the orientation of the B frame with 
respect to the GD frame. Due to the fact that four parameters are used to describe a rotation with three 
degrees of freedom, the constraint equation Cg + ef + 63 + 6§ = 1 must also be satisfied. 


2. Equations of Motion 

Newton’s laws of motion, expressed for a body moving relative to a planet-fixed reference frame, lead to the 
dynamic equation 


r + 2r2xf+r2xr2xf = 
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(30) 


as 


The quantity f may be expressed in the GC frame as f = —rezc- The velocity and acceleration follow 
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The angular velocity of the planet may be expanded as ^ — U {cosefexe ~ sin (j)ezc)- 

The force per unit mass acting on the vehicle, including the second zonal harmonic is given by 


F / 3/x J 2 . „ , \ ^ M 

— = 1 axe - sin 2<P ) bxc + aycByc + 1 aze + ^ “ 


(2 - 3cos^ e^c (33) 


m \ 2M I \ 

where axe, aye, and aze are the components of the geocentric representation of the applied acceleration. 
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where u, v, and w are the x, y, and z-axis components of the velocity in the GG frame, respectively, and 
axe, aye, aze are the geocentric representation of the accelerations acting at the body center of mass, not 
including the gravitational terms. 

The sensed accelerations at the vehicle center of gravity (denoted by axB, ays, and ozb) are transformed 
into the GC frame to be consistent with the equations of motion by means of 


(40) 
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The rotational kinematics of a vehicle relative to the GD frame may be written as 
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Note that the B frame accelerations at the IMU location axT, clyt, clzt and the B frame inertial angular 
rates, P, Q, and R are assumed to be known functions of time. The next section discusses the application 
of the filtering and smoothing methods to the Hyper-X Mach 10 test flight. 

III. Application to the Hyper-X Mach 10 Flight 

This section summarizes the results of applying the methods developed in the previous sections to the 
problem of estimating the trajectory of the Hyper-X vehicle on its third test flight, which occurred on 
November 16, 2004. This test flight was designed to reach a test condition of approximately Mach 9.6. The 
flight data preprocessing methods are discussed first. 

A. Flight Data Preprocessing 

This section summarizes the data preprocessing that were required in order to provide useful measurements 
as needed for the trajectory reconstruction process. The raw flight data were passed through a preprocessing 
tool, known as PREPR,^^ which consists of four main utilities. The first utility removes duplicate points in 
the data series which result from differences between the sensor data refresh rate and the telemetry sampling 
rate. The second utility allows users to interactively remove rogue points from the data series. Users select 
individual or groups of points to remove from a graph of the data series. The third utility samples the 
data series at a uniform time step, using a user-specified option of linear, cubic spline, or Hermite spline 
interpolation between available data points. The fourth utility gives the user the option to Alter the data 
series using Fourier series and optimal Weiner Altering method^^ available within the Systems Identification 
Program for Aircraft (SIDPAC).^^ A separate fifth utility is included which can be used to differentiate a 
data series once it has been passed through the previous four routines. The user is given the option to use 
backward difference or center divided difference methods. Alternatively, users can select to curve fit the data 
using either cubic or Hermite spline polynomials, and to individually differentiate each piecewise segment. 

1. Inertial Data 

Measurements of acceleration and angular rates were provided by QA2000 accelerometers and GG1320 Digital 
Laser Gyroscopes (DLG) contained within the Honeywell H-764G system, a unit with a background of 
numerous civilian and military applications. These measurement data were obtained from the on-board 
telemetry stream at a rate of 100 Hz and were first preprocessed by removing points which contained 
obviously bad or rogue data. Next, the data were edited to remove points which could be considered as 
repeated, stale, or non-updated frames. Gaps in the data due to telemetry dropouts or due to edited data 
were filled in by a piecewise cubic Hermite interpolating function, and were then Altered to remove the effects 
of high-frequency structural noise, measurement noise, and to smooth out the data from the effects of sensor 
quantization. Numerical derivatives of the inertial angular rates, as required for evaluation of the equations 
of motion, were calculated using a second-order centered divided difference scheme, except for the beginning 
and end points, which used first-order forward and backward divided differences, respectively. 

An example of raw data and processed inertial flight data is shown in Fig. 6, which shows the pitch rate 
(R-frame y-axis angular velocity measurement) during boost. The raw signal shows the evidence of a known 
8 Hz structural mode that was eliminated in the course of data pre-processing. 

Values of the process noise standard deviations used in connection with the inertial measurement data in 
the EKF formulation were taken from the Honeywell specifications given in Ref. 25. The inertial measurement 
data were also corrected for misalignment between the IMU mounting and the vehicle B frame, measured 
before flight and given in Ref. 26. Uncertainties in the preflight mounting measurements were provided and 
were taken into consideration when computing the uncertainties in the trajectory estimate. 

2. GPS Data 

GPS measurements were obtained as outputs from the H-764G GPS Embedded Module (GEM), which is a 
5-channel receiver, operating in Precision Positioning Service (PPS) mode at a rate of 1 Hz. The GPS data 
required little preprocessing, aside from extracting the data from the telemetry stream. There were several 
issues with the data itself, however. First, the east velocity channel saturated at approximately 27 sec. after 
the HXLV was dropped from the B-52 carrier vehicle, due to inadequate telemetry system bit assignments. 
The east velocity data unsaturated during the vehicle descent, at a time of approximately 688 sec. after 
drop. This saturation was a known problem and is due only to telemetry and is not a problem with the 
H-764G system. 
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Figure 6. Comparison of Raw and Processed Pitch Rate During Boost 


A second concern is that the GPS antenna carried on the HXRV itself was activated at the instant of the 
drop (previous GPS points from the HXRV telemetry stream were provided by the vehicle flight management 
computer, but were derived from a GPS antenna carried on the B-52), and as a result, suspected blockage 
or multi path effects may have corrupted the first several data points recorded after the drop. Therefore, 
the first several points were deleted from the data set and were not considered in the reconstruction process. 
Also, the GPS receiver lost lock on all satellites toward the end of the descent trajectory at a time of 
approximately 755 seconds after drop. The lock was regained at 780 seconds. During this period there were 
no substantial roll oscillations and the cause for the loss of lock is currently unknown. The data points in 
the telemetry stream during this time interval were removed from consideration for trajectory reconstruction 
purposes. 

B. Atmosphere Reconstruction 

The atmosphere reconstruction for the third flight of the Hyper-X vehicle is based on weather balloon data 
acquired both before and after the actual test flight. The atmosphere reconstruction is described in Ref. 27. 
The anchor point for the atmosphere reference was = 33.45°, 9ref = —121.0667°. Weather data at a 

point other than the location of the reference point must be transformed by making use of pressure gradient 
information to map the pressure to such a location. The temperature table is assumed to be invariant with 
respect to latitude and longitude, and the density follows directly from the application of the equation of 
state. 

The vehicle ground track relative to the atmosphere reference point is shown in Fig. 7. The drop, 
separation, and loss of signal points are also shown for reference. Glearly the evaluation of atmospheric 
data along the trajectory requires a large extrapolation from the region in which the data were originally 
captured, which represents a significant source of uncertainty. 

C. Results of the Reconstruction 

The results of the trajectory reconstruction for the Mach 10 flight of the Hyper-X vehicle is discussed in this 
section. The comparison of RMS errors between the reconstructed solution and the inertial solution with 
the GPS measurements is shown in Table 2. In all cases the trajectory reconstruction is a better fit for the 
GPS data than the pure inertial solution, as would be expected. 

An example of the improvement in the trajectory estimate gained by using the Kalman filter reconstruc- 
tion method can be seen in Fig. 8. This figure shows the north velocity solution during the boost phase, 
in which the reconstructed solution is shown in the solid curve, the on-board inertial solution is shown in 
the dashed curve, and the individual GPS solutions are shown with the square points. The drift the in the 
inertial solution is apparent, but is seen clearly in the residual plot, which shows the north velocity residuals 
during boost. The inertial solution residual shows a clear trend toward a value of approximately 1 ft /sec 
while the residuals of the reconstructed solution show no trend or bias. The maximum residual of the inertial 
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Figure 7. Ground Track Relative to Atmosphere Reference 


Table 2. Summary of RMS Errors 


Measurement 

EKF 

NAV 

Ratio 

Observations 

hgd^ ft 

33.52 

35.51 

1.059 

710 

(pgd, deg 

2.415-10-5 

7.890-10-4 

32.67 

710 

e, deg 

3.006-10-4 

1.441-10-5 

4.794 

710 

Ugd, ft /sec 

0.3979 

0.6425 

1.615 

710 

Vgd, ft/sec 

0.8183 

1.300 

1.589 

104 

Wgd, ft/sec 

0.7546 

0.8374 

1.110 

710 


solution is small compared with the magnitude of the total velocity, but nevertheless it is a systematic error 
that has been compensated for by using the filter/smoother approach presented in this paper. 

The next several subsections discuss the trajectory reconstruction and some comparisons with the on- 
board pure inertial solution during the boost, separation, and descent phases of the flight. Comparison 
against the preflight nominal trajectory and Monte-Carlo error bands are also shown. 

1. Boost Trajectory 

The altitude, Mach number, dynamic pressure, and angle of attack time histories during the boost phase 
of the trajectory are shown in Fig. 9. Here, the preflight Monte-Carlo average trajectory is shown in the 
dashed-dot line type, the Monte-Carlo +/- 3 standard deviations away from the nominal trajectory are 
shown in the dashed line type, the best-estimated trajectory is shown in the solid line type, and the +/- 
3 standard deviations of knowledge uncertainty are shown in the dotted line type. The pre-flight boost 
Monte-Carlo simulation is discussed in Ref. 4. 

Clearly these data indicate that the altitude and Mach profiles are well within the pre-flight predictions. 
The dynamic pressure is initially below the 3 ct pre-flight prediction, but the difference between the BET and 
the Monte-Carlo average rapidly decreases and for the remainder of the boost trajectory the dynamic pres- 
sure is essentially nominal. The initial deviation between the nominal dynamic pressure and reconstructed 
dynamic pressure is due to a difference in the range-reference atmosphere and the reconstructed atmosphere 
as well as a difference in the nominal and reconstructed velocity. 

The angle of attack time history starts near nominal but deviates briefly beyond the +/- 3 Monte-Carlo 
standard deviations in the transonic region and in the initial drop to the negative angle of attack region. 
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Figure 8. North Velocity During Boost 
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Figure 9. Trajectory Data During Boost 
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2. Separation 


The conditions at the separation between the HXLV and the HXRV are summarized in Table 3. The time 
of the commanded separation is 88.35 sec after the sensed drop time. The table lists the target separation 
condition and the range of acceptable deviation as specified in Ref. 28, along with the separation condition 
as indicated by the inertial navigation system (INS) with the day-of-flight atmosphere reconstruction. The 
last two columns give the separation condition determined by the trajectory reconstruction (BET) and the 
associated 1 cr uncertainties. 


Table 3. Flight 3 Separation Conditions 


Variable 

Target Requirements 

INS 

BET 

BET Uncertainty 

Time Since Drop, sec 

104.0 (Max) 

88.35 

88.35 

0.0 

Altitude, ft 

109,580 

109,601 

109,653 

6.4 

Mach number 

9.6 ± 0.6 

9.703 

9.702 

0.057 

Dynamic Pressure, Ib/ft^ 

1000 ± 210.0 

957 

954 

9.9 

Flight Path Angle, deg 

1.5 ± 1.0 

1.62 

1.66 

0.08 

Angle of Attack, deg 

0.0 ± 1.5 

-0.028 

-0.011 

0.043 

Angle of Sideslip, deg 

0.0 ± 1.0 

-0.20 

-0.20 

0.054 

Roll Angle, deg 

0.0 ± 15.0 

1.72 

1.67 

0.16 

Pitch Rate, deg/sec 

-0.18 ± 3.0 

0.0027 

0.414 

0.005 

Yaw Rate, deg/sec 

0.0 ± 3.0 

0.304 

0.292 

0.005 

Roll Rate, deg/sec 

0.0 ± 6.0 

0.032 

-0.05 

0.005 


Trajectory data during the separation of the vehicles is shown in Fig. 10. Similar to the boost trajectory 
data shown earlier, the pre-flight Monte-Carlo average trajectory is shown in the dashed-dot line type, the 
Monte-Carlo -f/- 3 standard deviations away from the nominal trajectory are shown in the dashed line 
type, the best-estimated trajectory is shown in the solid line type, and the -f/- 3 standard deviations of 
knowledge uncertainty are shown in the dotted line type. The Monte-Carlo separation simulation is discussed 
in Ref. 29 and background information on the stage separation system with flight data evaluation is provided 
in Ref. 30. All trajectory data are well within the pre-flight predictions, indicating an essentially perfect 
separation between the launch vehicle and the research vehicle. 

3. Engine Test Trajectory 

The scramjet engine test phase of the mission began with the opening of the engine cowl door at a time 
of approximately 2.5 sec. after the separation of the HXRV from the HXLV. The engine was fuelled for 
approximately 10 sec. The vehicle achieved a maximum Mach number during the scramjet engine test of 
approximately 9.6. The trajectory during this phase of the mission was again extremely close to the pre- 
flight predictions, but may not be shown due to sensitivity concerns. The scramjet engine test, including 
post-flight comparison with the pre- flight predictions, is described in Ref. 31. 

4- Descent Trajectory 

After completion of the scramjet engine test, the engine cowl door was closed and the vehicle then executed 
a pull-up maneuver to slow the descent. After completion of this maneuver, the vehicle was commanded to 
follow a descent trajectory that included a series of parameter identification maneuvers at each Mach number 
occurring from 8 to 2. The design and analysis of these parameter identification maneuvers is discussed in 
Refs. 32 and 33. 

The reconstructed trajectory data during the descent is shown in Fig. 11. Again, the pre-flight Monte- 
Carlo average trajectory is shown in the dashed-dot line type, the Monte-Carlo +/- standard deviations 
away from the nominal trajectory are shown in the dashed line type, the best-estimated trajectory is shown 
in the solid line type, and the -I-/- 3 standard deviations of knowledge uncertainty are shown in the dotted 
line type. The trajectory data are, for the most part, well within the pre-flight Monte-Carlo bounds. The 
main exception is the roll angle time history. Shortly after the closing of the engine cowl door the vehicle 
experienced some rocking motion. The cause of this motion is still under investigation.^^ 
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Figure 10. Trajectory Data During Separation 
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IV. Conclusions 


This paper has discussed the formulation and implementation of an Extended Kalman Filtering algorithm 
for use in reconstructing the flight path of the Hyper-X/X-43A vehicle, using in-flight inertial measurements 
in the dynamic model combined with redundant measurements of the vehicle state from independent sources, 
particularly from the Global Positioning System. The equations of motion were formulated in a local, non- 
inertial coordinate system, and included the effects of several systematic error sources. Equations relating 
the state to measured quantities were provided for a variety of data sources of interest to this problem. 
Data recorded during the Hyper-X Mach 10 flight on November 16, 2004, were first preprocessed and then 
processed using the extended Kalman filter, leading to an improvement over the on-board navigation solution. 
Results for the boost, separation, and descent phases of the flight were shown. 
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